Speeding up a single-molecule DNA device with a simple catalyst 
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Recently, several groups have designed and synthesized single-molecule devices based on DNA that 
can switch between different configurations in response to sequential addition of fuel DNA strands. 
There is considerable interest in improving the speed of these 'nanomotors'. One approach is the use 
of rationally designed DNA catalysts to promote hybridization of complementary oligonucleotides. 
A particularly simple and robust DNA device reported by Li and Tan is comprised of a single- 
strand 17-base oligomer that folds into a chair-like quadruplex structure. We have identified the 
key rate-limiting barrier in this device as the tendency for one of the fuel strands B to fold into the 
quadruplex configuration of the device strand. This seriously impedes the restoration reaction. We 
have designed a catalytic strand to inhibit the folding of B, and shown that the catalyst speeds up 
the restoration reaction by roughly a factor of 2. The catalyst remains effective even after repeated 
cycHng. 

PACS numbers: 87.15.He,81.07.Nb,81.16.Dn 



I. INTRODUCTION 

The operation of nanoscale devices based on designed 
DNA strands capable of switching between configura- 
tions in response to the addition of DNA fuel molecules 
in a cyclic manner has been demonstrated by several 
groups 0, 0, S S IE 0- The different configurations 
and the sequence recognition properties of DNA have 
been utilized to induce mechanical motion on a molecular 
scale. Mao et al. jl| exploited the transition between the 
B and Z forms of DNA molecules to induce atomic dis- 
placement of 20 to 60 A. The construction of nanomotors 
based on DNA molecules and branch migration has been 
pioneered by Yurke and his collaborators. These molecu- 
lar devices, dubbed "nanomotors" , exploit configuration 
changes induced by hybridization reactions. 

Yurke and collaborators have designed DNA tweez- 
ers 0, a DNA actuator Q, as well as a device that 
is switchable between three mechanical states 0. The 
transition between different states may be induced by a 
change in experimental conditions or by the addition of 
a "DNA fuel" that drives the cycle. Yan et al. @ have 
also reported a robust sequence-dependent rotary DNA 
device that flips between two well-defined stable config- 
urations. Progress in the design of single-molecule DNA 
devices with low molecular weight seems encouraging for 
future applications in nanotechnoglogy. A well-designed 
single-molecule DNA device should display a robust tran- 
sition between the straightened state and a folded state, 
examples of which are the G-quadruplexes 0, , and the 
hybridized hairpin structure 



II. CATALYTIC CONTROL OF THE DNA 
DEVICE 

Recently, Li and Tan ^ introduced a single-molecule 
device based on a short DNA oligomer. In their design, 
the 'motor' strand M is a 17-base single strand (5'-TGG 



TTG GTG TGG TTG GT-3') that folds up into a "chair" 
quadruplex structure in the presence of potassium ions 
(see Fig. QJi). The cycle is driven by the introduction 
of the 27-base fuel strand A (5'-GTA GTG CGG GAG 
GAA GGA GAG GAA GGA-3'), which consists of a 17- 
base segment that is the Watson-Grick complement of Af, 
and a 10-base overhang. Because the persistence length 
of double-stranded DNA is about 100 base pairs [iQj, M 
straightens out into a rigid double strand on hybridizing 
with A (Fig. ^p). It may then be restored to its initial 
chair configuration by the introduction of the restoration 
fuel B (5'-TGG TTG GTG TGG TTG GTG GGG GAG 
TAG-3'), which is the complement of A. The strand B 
first attaches itself to A via the 10-base overhang. Then 
it competes with M for binding with A through the pro- 
cess of branch migration ^|. When M falls off from A, 
it returns to the folded chair structure, completing the 
cycle. The complex AB (the waste product) does not 
participate further in the cycle. By successively adding 
A and then B, the motor strand M can be repeatedly cy- 
cled between its straightened and folded configurations. 

As shown in Fig. the strand M is labeled at 

its 5' and 3' ends with the fiuorophore 6-FAM (5'- 
fiuorescein phosphoramidite) and the quencher DAB- 
GYL [4-(4'-dimethylaminophenylazo) benzoic acid], re- 
spectively. At room temperature, the excitation and 
emission peak wavelength of the 6-FAM dye are 492 nm 
and 518 nm, respectively. The molecule DABGYL is a 
wide-spectrum quencher that quenches efficiently the flu- 
orescence of 6-FAM. Because fluorescence is unaffected 
(strongly quenched) when the ends are far apart (in close 
proximity) , the fluorescence intensity is a reliable indica- 
tor of the conflguration change in M . 

In the device of Li and Tan Q , the fuel strand B has a 
segment that is identical with M (by complementarity). 
Hence it shares the latter 's tendency to fold. Folding of 
a fraction of B^s explains why the restoration reaction is 
significantly slower than the straightening reaction. The 
method of mutated fuel was developed by introducing 
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mismatches between the motor strand and straightening 
fuel 1^, |9(| . This partiaUy solves the problem, but at the 
cost of slowing the straightening reaction. In a recent 
experiment, Turberfield et al. 12] used loop complexes 
to inhibit the hybridization of complementary oligonu- 
cleotides. They found that rationally designed DNA cat- 
alysts are effective in promoting the hybridization. 

Here we report a catalyst that specifically applies to 
the DNA device of Li and Tan Q. The catalyst, which 
works by binding to the fuel strand B in order to in- 
hibit its folding, consists of a short strand C (5'-CGC 
GAG GA-3') that hybridizes to 2 of the 8 guanines in 
B that form the quadruplex. In addition, C hybridizes 
with other nearby bases in the overhang of B to stabi- 
lize the binding. With 1 M monovalent cation, at the 
concentration of 1 /iM, the melting temperature of BC 
is 45.6°G 13], which is higher than room temperature 
at which the device is operated. (In the absence of the 
overhang C hybridizes with only 4 bases in M . As the 
melting temperature of the complex CM is below room 
temperature, C does not form a stable duplex with M.) 
If the folding of B is the reason for slow restoration, the 
catalyst C should accelerate the reaction (Fig. ^). 



III. EXPERIMENTAL RESULTS 

All oligonucleotides used in the experiments were syn- 
thesized and purified by Integrated DNA Technologies, 
Inc. (IDT). Stock solutions were prepared by resus- 
pending the lyophilized oligonucleotides in sodium phos- 
phate/sodium chloride (SPSG) buffer (pH 6.5, 50 mM 
Na2HP04, 1 M NaGl) at the concentration of 100 /iM. 
The experiments were performed at 1 /LtM concentra- 
tion at 20°G by diluting stock solution of M 100 fold 
in KSPSG (pH 6.5, 40 mM KGl, 50 mM Na2HP04, 1 
M NaGl) buffer. The stock solutions of strand A and B 
are added and mixed by rapid pipetting up and down for 
-10 s. 

Time-based fluorescence intensity curves were taken on 
a custom-designed spectrofluorometer (Photon Technol- 
ogy International), using a Hellma 105.254-QS cuvette. 
The excitation wavelength was set at 492 nm and the 
emission wavelength at 518 nm. The excitation light 
intensity was kept low so that during the experiment 
the fluorescence bleaching effect was negligible compared 
with the effect caused by configuration changes in M . 

Figure 121 compares the time dependence of the fluores- 
cence intensity with the catalyst present (Sample a) or 
absent (&) as the device is cycled by successive addition of 
fuel A and B. In both traces, the straightening reaction 
on addition of A is too rapid to be resolved. By con- 
trast, the restoration reaction following the addition of 
B is slower. Gomparing the two curves, we observe that 
the addition of C does not affect the steep rise when A 
is added (within our resolution capability), but signifi- 
cantly accelerates the restoration reaction following the 
addition of B. 



The vertical jumps in both traces when A or _B is added 
(at time t — and 120 s, respectively) is an artifact 
of our experimental procedure. Right before A ov B \s 
added, the shutter of the spectrofluorometer is turned 
off for a short 'blind' interval (Ai ~ 10 s) to protect 
the photo-sensors from strong light. During this interval, 
the solution is repeatedly pipetted to improve mixing, 
but data are not recorded. The relatively fast reaction 
time (~ 100 s) in the Li- Tan device makes this 'blind' 
interval much more noticeable than in motors with longer 
reaction times. However, the numerical fitting procedure 
(Appendix) allows us to interpolate the data within At. 
(During At, the highly inhomogeneous concentration of 
B prior to mixing seems to accelerate the reaction greatly 
in some regions of the solution. Our simulations favor a 
longer effective At than the 10 s when the spectrometer 
is actually in the off state.) 

The effect of varying the concentration of the catalyst 
C on the restoration reaction has been investigated in 
detail in an ensemble of 24 samples (at 8 different con- 
centrations). To minimize the experimental uncertain- 
ties, these samples were all prepared at the same time, 
kept wrapped in aluminum foil and stored at 4°G until 
the measurement. In Fig. three curves a, h and c 
at successively higher concentration of C are shown in 
linear scale. Glearly, the rate of decrease is accelerated 
with increasing concentration of C. In each case, rather 
close fits are obtained using the model described below 
(solid lines). The fits allow interpolation within the blind 
interval just after B is added. 

Although numerical integration of the rate equations is 
necessary to describe the curves accurately, we note that 
at early times (semilog plot in Fig. I^a), they are well- 
described by a single exponential term. The straight-line 
behavior in early times determines an effective time con- 
stant T which we have used to characterize the effect of 
increasing C (at latter times the experimental curves de- 
viate noticeably from simple exponential decay). The 
dependence of r on the concentration of C measured in 
the 24 samples is displayed in Fig. 0] which is our central 
experimental finding. At low concentrations, the catalyst 
speeds up the reaction rate dramatically but the improve- 
ment saturates at moderately high concentrations. [We 
emphasize that the absence of data at the very early in- 
terval At ('vertical jump' discussed above) does not pose 
a problem for the semi-log plot in Fig. as it merely 
shifts the origin of t. Using the semi-log plot, the ef- 
fective time constant r is directly inferred from the raw 
data. Hence the variation of r shown in Fig. 2|is inde- 
pendent of our model.] 

The catalyst remains efficient under repeated cycling. 
As shown in Fig. [S] the restoration reaction with catalyst 
is reproducibly faster than that without catalyst. 
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IV. ANALYSIS 

The operation of this DNA nano-dcvicc is driven by the 
hybridization of Watson- Crick complementary sections of 
DNA. In the absence of the catalyst, two processes are 
involved in the restoration reaction. First, B hybridizes 
with the single-stranded overhang section of A. Following 
the opening of B, random branch migration along the 
motor-strand section of MA occurs until, either the waste 
product AB is formed, and M returns to its folded state, 
or the branch point returns to its initial point and B 
dissociates from the complex (or branch migration starts 
anew) . 

As observed in our experiment, the restoration is al- 
ways much slower than the straightening. In branch mi- 
gration, the random walk process is associated with the 
temporary breaking of base pair bonds by thermal fluc- 
tuations. When uninhibited, a random walk over short 
segments (10-15 base pairs) is expected to take about 
15-20 ms 1131, which is much shorter than the observed 
reaction time. Hence something else is impeding restora- 
tion. An important factor here is the tendency of the 
unhybridized B to collapse into its folded state. We rea- 
son that the folding of B strongly inhibits the initiation 
of branch migration right after the overhang sections of 
B and A become hybridized (see Fig. ^Jd). Hence this is 
most likely the reason for the relatively long restoration 
time of the cycle. 

In the "all-or-none" model of Morrison and Stols |Q| , 
the hybridization reactions of DNA is described by the 
equation 

S + S% SS, (1) 

where S and S are a pair of complementary DNA strands, 
and fc+ and fc_ are the forward-reaction and off rate con- 
stants, respectively. 

Our model is then described by the following reactions 



B + C 


kg ^ 


BC 


(2) 


MA + B 


fci^ 


MAB 


(3) 


MAB 


k2^ 


M + AB 


(4) 


MA + BC 


k3^ 


MAB' + C 


(5) 


MAB' 


k4 


M + AB, 


(6) 



where fcg, • ■ ■ , A:4 are the corresponding forward reaction 
rate constants (the subscript + is suppressed hereafter). 
In the Appendix, we show that the corresponding off re- 
action rate constants are negligible. 

Equations I2I6I which describe 2 parallel processes for 
the production of the final products M and AB ('sinks') 
starting with the products MA and B ('sources'), may 
be represented by the simple diagram in Fig. El The 
'slow' process (Eqs. |21 and ^ is independent of C, and 
involves the hybridization of MA and B to form MAB 



which subsequently decays to M and AB (upper branch 
in Fig. EJ. Its rate constants, ki and fe, are readily 
obtained by fitting the measured decay curves in which 
Co = 0. (see Appendix). 

The 'fast' process, which is catalyzed by C, is initi- 
ated by the formation of the duplex BC (Eqs. [3 , El 
and El lower branch in Fig. EJ. Hybridization of BC 
and AIA produces the unstable 4-strand complex MABC 
which quickly decays to MAB' and C, and finally to 
M and AB. For clarity, we have not shown MABC in 
Fig. El We distinguish MAB' from MAB (in the up- 
per branch) because in our model branch migration is 
assumed to occur more readily in the former, so it has a 
faster decay rate than MAB. The observed fluorescence 
intensity is proportional to the sum of all non-free M, 
[MA] + [MAB] + [MAB']. 

As implied in Fig. El the fast and slow reactions com- 
pete for B. Our experiment shows that the competition 
is strongly tilted towards the lower branch as soon as the 
catalyst concentration cq is increased from zero (steep 
drop in r in Fig. 0}. The saturation of t at large cq 
also implies that the slow process becomes insignificant 
for Co > 0.5 fiM. Hence the curves at 1.5 fxM are deter- 
mined by fco, fcs and /c4. It turns out that /c4 is the most 
important parameter determining the observed decay in 
the high-Co regime. The "bottleneck" at the node MAB' 
resulting from the relatively small value of compared 
with k^ makes the decay rate rather insensitive to the 
latter (the lines in Fig. El are drawn with widths roughly 
proportional to the rate constants). By fitting the curves 
in the 2 extreme limits, co = and co = 1.5 /iM (Ap- 
pendix), we find the values of /ci, fc2 and to be 

fci = (0.21 ±0.04) fiM~h-\ (7) 
k2 = (8.9 ± 1.0) X 10"^ s-\ (8) 
fc4 = (3.4 ±0.1) X IQ-^ (9) 

[We remark that the second-order rate constants fco, fci 
and fca (reactions on the left half of Fig. EJ carry an ex- 
tra factor of concentration compared with the first-order 
constants fc2 and fc4 (right half). Since the concentra- 
tions used here are nominally fiM, we have expressed the 
former in units fj,M~^ s~^, and the latter in s~^. In these 
units, the numbers quoted provide a more meaningful 
comparison of the reaction rates on the left half with 
those on the right.] 

While the remaining rate constants fco a-nd k^ only 
weakly affect the decay rates of the large-co curves, they 
obviously infiucnce the strong competition for B when co 
is small. Hence the observed curve of t vs. co in Fig. ^ 
may be used to estimate 

fco = (2±1) ^lM-h-\ (10) 

fca = (2.5 ±0.5) ^M-^s-^ (11) 

Although the uncertainties in fco and fca are relatively 
large, our fits exclude values of fco < 1 /iM~^s~^ and 
fcs < 2 ^M~^s~^. The large value of fcs relative to fc4 
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leads to the bottleneck at the node MAB' mentioned 
above. 

The ratio of the reaction rate constants ki/k2 ~ 3.8 
implies that the restoring reaction Eq. ^is significantly 
improved with the addition of C. The large value of 
(~ lOfci) implies that C is highly efficient in catalyzing 
the formation of the complex MAB' compared with the 
formation of MAB in the catalyst-free solution. This 
presumably reflects the inhibition of the folding of B. 



V. SUMMARY 

In summary, we have designed a DNA catalyst strand 
that significantly accelerates the restoration reaction 
time in the device proposed by Li and Tan. Recently, 
Dittmer, Renter and Simmel |16| reported an interest- 
ing application of a DNA aptamer similar to this de- 
vice. They successfully used the machine to precisely 
control the concentration of thrombin protein in solution 
between a depleted and an enriched state. Our catalyst 
may be useful in their application. A model incorporat- 
ing the reaction rate constants of the intermediate prod- 
uct helps us to understand the rate-limiting features in 
the operation of this single-molecule device. As noted 
in Ref. the catalyst itself is also a DNA machine 

that obtains energy by catalyzing the restoration reac- 
tion. Such machines can run until the supply of unre- 
acted fuel is exhausted. 



Appendix 

We discuss in more detail our analysis, and the opti- 
mization procedure. First, we justify the neglect in Eqs. 
121 and 01 of the off reaction-rate constants fc_ (defined in 
Eq. HJ. We have 



K 



X IM 



^-AG/RT 



(12) 



where AG is the reaction Gibbs free energy, R the molar 
gas constant and T the absolute temperature. The Gibbs 
free energies are estimated from the data of Peyret and 
SantaLucia In a solution at 20°C with 1 M mono- 
valent cation, the Gibbs free energies are AG ma = - 
27.55 kcal/mol (for hybridization of M and A), AGab' 
= -16.13 kcal/mol (of A and B) and AGbc = -13.40 
kcal/mol {B and C). In our experiment the nominal val- 
ues of [S]o are 1 /iM, which implies that fc+[5]o/fc_ ^ 1 
in each reaction (the smallest value - fc+[5]o/fc_ = 9,900 
- is obtained in the reaction Eq. {21 . Hence, in compar- 
ison with the forward rate constants, we may ignore fc_ 
in Eqs. 1^ and 151 [l5l| . In the reactions Eqs. I4I6I the re- 
verse reactions, which involve the dehybridization of the 
stable complexes MAB' and AB, have a negligibly small 
probability as well. 



In terms of the the forward rate constants, the time 
derivatives of the populations are 



d[BC] 
dt 

d[MAB] 
di 

d[MAB'] 
Jt 

d[M] 
dt 



= ko[B][C] (13) 

^ ki[MA][B]- k2[MAB] (14) 

^ k3[MA][BC]-k4MAB'] (15) 

^ k2[MAB] + k4MAB']. (16) 



In the experiment, the total amount of Af , A and B 
add up to 1.0 ^M, while the amount of C is cq. Hence 
by mass conservation, we have (with the assumption [A] 
= 0) 

[M] + [MA] + [MAB] + [MAB'] = 1.0 /iM (17) 



[MA] + [MAB] + [MAB'] + [AB] = 1.0 fiM (18) 



[B] + [BC] + [MAB] + [MAB'] + [AB] = 1.0 fiM (19) 



[C] + [BC] = Co. 



(20) 



Equations ini and UHl imply [M] = [AB]. 

For a given co, we have as input the measured fluores- 
cence intensity curve (normalized to 1 at t = before the 
addition of B) 



fit) 



Ijt) - I. 
1(0) - lo 



(21) 



(with the caveat that data is missing in the blind interval 
At). For a set of starting values of ki's we integrate 
numerically the rate equations to determine the curve of 
1 — [M]{t) (hereafter, 1 means 1.0 fiM). A least-squares 
fit is used to fit the curve of 1 — [M] to f{t) as follows. 

We note that the rate equations Eqs. I14ll6l simplifv in 
the 2 limits of zero concentration and high concentration 
of C. In the first case, we have [BC] = [MAB'] = 0, 
which gives (from Eqs. 1171 and lT^ 

[B] = [MA] = 1-[M]-[MAB] (cq 0). (22) 

Eg nations 1 1 41 1 61 reduce to 

d[MAB] 



dt 

d[M] 
dt 



= ki{l - [M] - [MAB]f - k2[MAB]{23) 
= k2[MAB]. (24) 



In the curve-fitting procedure, we assume seed values 
for ki and fc2, and integrate Eqs. 1231 and OH numericallv 
to obtain the curve of y{t) = 1 — [M](f), which is then 
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compared with the measured f{tobs)- Here, tobs = t — At 
where At is the shght shift in the origin of the time 
axis caused by the mixing step, as explained in Sec. IIIII 
The values of ki, k2 and At are successively refined un- 
til convergence is attained. The convergence is faster 
if we regard y and / as the independent variables, and 
perform the least-squares fit on the inverted curves t(?/) 
and tobsif) + At. For the 3 samples with cq = 0, the 
optimization process yields ki = (0.21 ± 0.4) /iM~^s~^, 
k2 = (8.9 ± 1.0) X 10-3 s-i and A<o = 56 s. 

By fitting to the curves with cq = 1.5 /iM, we find 
ki = (3.4 ± 0.1) X 10-2 s-\ and Ati = 21 s. (As noted, 
these large-Co curves are highly sensitive to variations 



in ^4, but only weakly sensitive to variations in ko and 
fca.) With the values of fci, ^2 and k^ determined, we 
can calculate the effective time constant r vs. cq using 
assumed values of /cq and ^3. Finally, variation of these 
values allows us to optimize the fit of r vs. cq to the 
measurements in Fig. 0] The best values are reported in 
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(a) 



/ DABCYL 
6-FAM 





FIG. 1: Operation of the single-molecule DNA device, 
(a) The motor strand M in solution. M is labeled at 
end 5' with the fluorophore 6-FAM (5'-fiuorescein phospho- 
ramidite) and at end 3' with the quencher DABCYL [4-(4'- 
dimethylaminophenylazo) benzoic acid]. When folded into 
the chair structure, close proximity of the two ends quenches 
fluorescence. Dashed lines indicate the two G-quadruplexes in 
folded configuration, (b) Reaction cycle of device introduced 
by Li and Tan 0. In the 'straightened' configuration of MA, 
the fluorescence intensity is high. The intensity decreases as 
folding is gradually restored with introduction of the restora- 
tion fuel B. Undesired folding of B slows the restoration reac- 
tion, (c) Reaction cycle with added catalyst. The catalyst C 
is designed to bind to B to inhibit folding. It also helps to ini- 
tiate the branch migration and accelerate the restoration re- 
action. For substoichiometric cq both (b) and (c) occur when 
B is added. For stoichiometric and above-stoichiometric co 
(c) dominates the restoration reaction. 
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FIG. 2: A reaction cycle of the single DNA molecule motors 

with (Sample a) and without (Sample 6) the catalyst. High 
fluorescence signal corresponds to having the ends of M far 
ajjart while low fluorescence means that the ends are in close 
proximity (M folded). In each curve, the sample is 100 /A 
of the motor strand M at the concentration 1 fjM in KSPSC 
buffer at t = 0. At time t = 30 s, 1 /xl of 100 /uM C stock 
solution was added to Sample a, while 1 fjl of SPSC buffer was 
added to the control sample b (i.e. the catalyst concentration 
is ~ 1.0 ijM in a and in b). Note that the catalyst accelerates 
the restoration reaction without affecting the straightening 
reaction. Right before or B is added, the shutter of the 
spectrofluorometer was closed (for ~10 s) and the solution 
rapidly pipetted. This accounts for the discontinuity at t = 
60 and 120 s (see text). 




time (s) 



FIG. 3: Time-dependence of 6-FAM fluorescence during 
the restoration reaction. The background quenched inten- 
sity from folded M has been subtracted, so the fluores- 
cence intensity is proportional to the sum of all non-free M, 
[MA] + [MAB] + [MAB']. The intensity before addition of 
B has been normalized to 1. The solid lines are fits to the 
Eqs. I2I6I (see Appendix) . The calculated curves with opti- 
mized reaction rate constants fci ■ • • ^4 interpolate within the 
'blind' intervals At. The catalyst concentration in Samples a, 
b and c are 0, 0.1 and 1.5 fiM C, respectively. In the inset, the 
same data are displayed in semi- log scale. The initial straight 
portion is used to determine an effective time constant r. 
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FIG. 4: Restoration reaction time constant as a function of 
catalyst concentration ([M] = [A] = [B] = 1.0 iiM). As noted 
in the text, values of t shown here are inferred directly from 
semi-log plots of the raw data, and are independent of fits to 
the model. The best fit to t vs. co (solid curve) is used to 
determine the reaction rate constant ks (Appendix). 




FIG. 5: Cycling the DNA device with 1.0 /jM of catalyst 
C (solid curve) and without catalyst (dotted curve). In each 
cycle, high (low) fluorescence signal indicates straightening 
(folding) of the motor strands M. The catalyst remains effi- 
cient for several cycles. 



FIG. 6: Schematic drawing of the reaction paths in Eqs. 12161 
The duplex AIA, fuel strand B and catalyst C are 'sources' 
on the left while the motor strand M and waste product AB 
are sinks on the right. In the absence of C, the reaction 
proceeds along the upper branch (with rate constants ki and 
^2). Introducing C turns on the process on the lower branch. 
For small co(< 0.5 fiM), the 2 reactions strongly compete for 
B, whereas for cq > 0.5 ^M, the lower branch dominates. 
The thickness of the connecting lines nominally reflects the 
rate constants. In both branches, the small values of the rate 
constants k2 and ^4 lead to bottle-necks in the intermediate 
products MAB and MAB'. 



